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We generalize results of lattice QCD to determine the spin-dependent symmetries and factorization 
properties of meson production in OZI allowed processes. This explains some conjectures previously 
made in the literature about axial meson decays and gives predictions for exclusive decays of vector 
charmonia, including ways of establishing the structure of F(4260) and F(4325) from their S- 
wave decays. Factorization gives a selection rule which forbids e + e~ — ► D* Di near threshold with 
the tensor meson in helicity 2. The relations among amplitudes for double charmonia production 
e + e~ — > i/>Xo,l.2 are expected to differ from the analagous relations among light flavour production 
OO ' such as e + e~ — > Lofa \ 2. 
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I. INTRODUCTION 



The dynamics of strong decay amplitudes are poorly understood. Definitive answers are not known to questions as 
basic as: (i) are the qq created in an OZI allowed decay spin singlet or spin triplet; (ii) what is their overall J PC ; 
^ (— I (iii) are the qq created from the energy of the strong confinement field, or from a hard gluon? It is our purpose 
O i 1 in this paper to address these questions. We shall show that results from lattice QCD imply that light qq pair has 
' I spin 1 with an effective factorization of constituent spin and orbital degrees of freedom such that the qq pair in the 
initial meson are passive spectators. By contrast, if heavy flavours are created, as in e + e~ — * ip + x,j, factorization is 
I broken with spin and momentum transferred from the initial cc to the created pair, such as by a hard gluon. This 
implies a radically different spin dependence of amplitudes, and of angular distributions, in analogue processes such 
as e + e~ — * cj/2 relative to e + e~ — ► ipX2- 

While Lattice QCD is now a mature guide for the masses of glucballs and hybrids, at least in the quenched 
approximation [l|, |2j, it is not yet mature enough to determine hadronic decays extensively. Consequently, at a 
fundamental level the dynamics of such decays are not yet established. Flux tube models of both spectra [!, H[ and 
decays have been developed, in part stimulated by attempts to model the lattice, and lattice work has confirmed 
I ' their spectroscopy 0, 3 ■ The lattice is now beginning to confirm aspects of the flux tube model for some decays: 
. specifically the lattice QCD studies of the decays of hybrid 1 h — » 7r6i and 7r/i Q show quantitative features that were 
anticipated in flux-tube models 0, Q , and in ref 0] we showed that these approaches exhibit remarkable agreement 
t^- , when compared under the same kinematic conditions. Specifically, for S-wave decay amplitudes at zero-recoil the 
■ results are consistent with 



0(7T1 -> 7T + fclpPl]) = 2a(7T! -» 7T + /l[ 3 Pl]). (1) 

where 7Ti denotes the first gluonic excited hybrid with J p = 1 ^. 

In section II we describe the underlying assumptions of the factorization hypothesis: (i) the hadrons' spins, j, 
separate into two parts - the intrinsic spins of the constituents s and a residual component that transforms as angular 
momentum I, (ii) the I and s degrees of freedom act independently throughout the transition ("factorization"), (iii) 
the qq pair produced has spin 1. In section IIA we demonstrate that the above ratio is immediate within the 
factorization hypothesis. We identify further implications of factorization for the decays of axial and vector mesons in 
sections IIB and IIC; the former confirms and explains a conjecture of [i| and the latter has implications for charmed 
meson production from a -0( 3 Si) initial state. A helicity selection rule is derived implying that in e + e~ — > D* D2 near 
threshold or from a ?/>( 3 Si) initial state the D2 cannot be produced with helicity two. In section IID we apply the 
results of factorization to the decays of cc to charmed mesons near threshold in relative S'-waves, and identify ways 
to distinguish between hybrid and conventional interpretations of enigmatic ip-like states such as F(4260) or F(4325). 
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In section III wc discuss application to e + e~ annihilation for light and heavy flavours. We show that decays triggered 
by hard gluons violate factorization, and that emerging data on e + e~ — > ip + Xj appear to support this. We propose 
tests for factorization and hard gluon production mechanisms in e + e~ — > i\) + Xj near threshold. 

II. FACTORIZATION: FORMULATION AND APPLICATION 

Our approach to strong decays is similar in spirit to what was done in past decades for electromagnetic and 
current induced transitions among hadrons, known variously as Melosh transformation or more generally "single 
quark transition algebra" [l^l ■ The empirically successful hypothesis there was that the interaction of a current with 
a single quark triggers a transition, all other constituents being passive spectators. That led to algebraic relations 
among amplitudes, which arose from the Clebsch-Gordan coefficients coupling the orbital, spin and total angular 
momentum projections l z , s z to j z for the initial and final hadrons, and the l z , s z algebraic transformation properties 
of the transition operators. While the relative strengths of the reduced matrix elements associated with each transition 
operator are in this general approach undetermined, the experimentally accessible range of hclicity amplitudes for 
photo-excitation of proton and neutron targets to different resonances within a supcrmultiplct led to experimentally 
testable relations among various amplitudes Within the hypothesis of I, s factorization, analogous relations arise 
for strong decays. Specific models have implicitly assumed such factorization 0, [|| [l2|, EH, EII; we shall see that the 
results from lattice QCD suggest that this property is realised in decays at least for light flavours. 

We consider the decay process of mesons M — » Mi + M 2 . In meson M with spin j, the qq have spin s and residual 
angular momentum I. We illustrate the structure of the amplitude for the particular case of a flux-tube that breaks 
to form a new qq in a 3 Po configuration leading to a pair of mesons Mi, M2 with spins Ji,j2 respectively, and their 
qq having s\,h, S2,h- The generalisation will be immediate. 

The width for the decay of a meson into a pair of mesons involves a sum over couplings of j\ ® 32 to j'12 and relative 
partial waves L: 

T(slj -> Sihh + s 2 hh) ~ ^2((((si ® Ji)* ® (s 2 8> h)j 2 )h 2 ® L)j\\a ■ V||(s ® l)j) 2 (2) 

where a transforms as a vector in spin-space and V acts on the spatial (orbital and radial) degrees of freedom. 
Usually at th is p oint specific wavefunctions are assumed and nonrelativistic expressions calculated for the amplitudes 
[S HI, Ell [HI, Q3| ■ However, this introduces model dependence and obscures the more general underlying properties. 
Instead we shall factor the amplitude in such a way that the spin and space parts are separated [la ], expressing all 
decay amplitudes as linear combinations of model-dependent spatial amplitudes, which in the present work are left 
general. 

The first step is to separate the spin and space degrees of freedom of the final state. The final state bra 

((((si ® h) n ® (s 2 ® h) n )j 12 ® L)j\ (3) 

is recoupled to form states of good (s\2,l\2,lf)'- 

(((si <8> s 2 ) Sl2 ® ((h ® h)i 12 ® L)i f )j\ (4) 

which involves a product of 6-j and 9-j coefficients. The spin and space parts then factorise and can be isolated. The 
result is 

((((*i ® h)n ® (52 ® h) n ) n2 ® L)j\\a ■ V\\(s ® l) 3 ) = ^ {-) s+L+Sl ' +h2+ll ^i s s l2 i l2nn3l2 

si 2 h 2 lf 



si h ji 

S2 h j 
Sl2 ll2 jl2 



S2 h 32 \\ Sl j> l f t/jlT i f )\(( s i ^^iHl^^i^^^i^HVIIO (5) 



with 



n Qb ... = y/(2a+l)(2b+l)... (6) 



3 



The spin part is a 9-j coefficient along with appropriate counting factors 

((*i®* 2 ).JMI*> = (-) s+Si n lsSlS2Sl J | \ s 2 \. (7) 

[ S 1 S12 J 

The 9-j coefficient in cqn. (J7J is zero for s± = s 2 = s = 0; this is the well known spin singlet selection rule and is a 
consequence of the orthogonality of the spin wavefunctions. Note in the above a phase of (—1) has been included for 
the permutation of quark and antiquark operators [j| , and the expression ([5]) is equivalent to that in Ref. . 

The assumption driving the expansion of eqn. ([5]) is that the angular momentum and spin quantum numbers 
factorize and that the decay operator is overall scalar with a spin triplet part. The angular momentum algebra makes 
no reference to the spatial part of the operator, hence the linear combinations for 3 Po and 3 Si decay models, driven 
by operators a ■ V and a ■ f respectively, are the same. In a 3 Si model the spatial contraction involves f - the unit 
vector in the relative coordinate of the initial meson's qq . In constituent gluon models the spin-dependence of 
the qq creation again is via a a operator, while the spatial contraction depends on the specific model wavefunctions. 
In general, for any specific model there will be differing spatial dependence but the overall spin coupling coefficients 
are identical. 

Eqn. ([5]) expresses full decay amplitudes Oj 12 l as linear combinations of model dependent spatial amplitudes Ai 12 n f 
of the form 

Ai 12 u f = (((h®l2) h2 ®L) lf \M\l) (8) 

The expansion applies to all partial waves L allowed by the conservation of angular momentum, including those 
which are parity-forbidden for a given set of spatial quantum numbers. It is the spatial matrix element itself which 
ensures the conservation of parity; this is verified in the expressions of Ref. [l5| where the spatial matrix elements 
for the production and decay of conventional and hybrid mesons are presented. Thus, for instance, decays of the type 
3 Pi — Pi 1 So are allowed in S-, P- and D-wave in general; if the initial 3 Pi is a conventional 1 ++ or hybrid l" 1 the 
S- and D-wave amplitudes vanish, whereas if the initial 3 Pi is a hybrid 1 ^ the P-wave amplitude vanishes. 

The approach taken here, however, is to exploit the relationships between the decay amplitudes a J12 £ leaving the 
spatial amplitudes Ai 12 u f undetermined. These spatial amplitudes depend on the decay momentum and the spatial 
wavefunctions; thus, in the limit that the spatial wavefunctions of the mesons under comparison are the same and the 
momenta are the same, the expansion of Eqn. ([5]) relates decay amplitudes among families of states sharing the same 
spatial quantum numbers but having different spin and total angular momentum. If for a given partial wave L there 
is only one spatial matrix element of the form ([5]), which we denote Al, there are direct relations among amplitudes 
for states with different angular momentum quantum numbers. Three such cases are immediate: 

1. one of the final states has orbital angular momentum zero (1 2 = 0) and the decay is in relative S'-wave [L = 0); 
thus If = I12 = h and the amplitude is expressed in terms of a single matrix element As; 

2. both final states have orbital angular momentum zero (l\ = I2 = 0); thus I12 = and If = L and the amplitude 
in a partial wave L can be expressed in terms of a single matrix element Al,; 

3. the initial state and one of the final states have orbital angular momentum zero {I = I2 = 0), thus I12 = h and 
If = 1 and the amplitude in a partial wave L can be expressed in terms of a single matrix clement Al . 

We now examine each of these three cases in turn with specific examples. In section Hi Al the S'-wave hybrid decays 
7Ti — > &i7T and 7Ti — * /i7r are shown to match results from lattice QCD and thereby to reveal significant information 
about the underlying dynamics (case 1). In section Til B[ a\ — > pit and b\ — > lot: arc examples of case 2; the analysis 
verifies a conjecture that was made elsewhere [j| and establishes its origin. In section Til CI case 3 is applied to derive 
relations among decays of the type 3 Si — > ( 1 Pi; 3 Pj) + ( 1 So; 3 Si). A new selection rule is derived and the possibility 
of testing it in the context of e + e~ annihilation producing flavoured and flavourless states is discussed. In section IID 
we discuss ways of using these results to distinguish hybrid cc from 3 Si or 3 Di ip states. 

A. S'-wave decays of hybrid meson 7Ti 

An immediate example of this factorization is the S-wave decays of the hybrid meson -k\ — > bin or fin. The tt\ has 
I ^ quantum numbers and j = s = I = 10, Q. In flux tube models and non relativistic constituent gluon models (l6j 
the I = 1 is explicit; in cavity and bag models it is implicit in the definition of the TE gluon mode which transforms 
as f x e (e.g. see eq. (2.22) in [l7| and applied to hybrid decays in [HI)- The final states have l\ = ji = 1 and 
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S2 = 1-2 = 32 = 0, differing in the spin quantum number s\ = (6i) and s\ = 1 (/i). For decays in S-wave there is 
only one matrix element of the form l[8]). having l 12 = If = 1. In the summation over l 12 and If, this constraint is 
enforced by zeroes in the 9- and 6-j coefficients which reduce to delta functions. This reduces the expansion of cqn. 
([5]) to a simpler form and the amplitude for the initial state with spin s± is given by 



-(■r^-P.+'fW-Xfii}^ V 



1/2 1/2 S! 

1/2 



x A s (9) 



where here As is the spatial matrix element. Thus for the 1 Pi and 3 Pi modes 

A s 



as( 3 Pi - 1 Pi ^o) = --| (10) 
a s ( 3 Pl - 3 Pl %) = (11) 

The flavour overlaps for 7Ti — * b\i\ and fi(nn)n are \/2/3 and yT/3, so that the ratio of amplitudes is 

as(7Ti -> &i7r) 



as(i"i ->• /itt) 



= 2 (12) 



which underwrites the result eq.(l) as found also in lattice QCD. The essential feature here is the factorization of 
spin and space, and the assumption that the created qq are spin-triplet. Note that the qq creation with quark-spin 
1 now appears explicitly in 9-j and 6-j symbols; if the created pair has spin the final 9-j symbol for the irbi mode, 
equation |7|). has a zero in the bottom left corner, and since 



I I ) 
1 



the decay m — > tt + &i( 1 Pi) would vanish. If spin is conserved, an initial state with S = 1 can only decay to a pair of 
S = states if S — 1 qq is present, hence the need for pair-creation to be spin-triplet for a non-zero amplitude. 

Thus the results of lattice QCD, at least when applied to the decays of a hybrid mesonQ, follow if the amplitude 
factors in space and spin, with the qq pair creation being spin-triplet and an overall scalar. This does not distinguish 
3 Po from 3 Si decay models. 



B. S and D wave decays of axial mesons 

Ackleh et al@ noted that the ratio of the D / S'-wave amplitude ratios for b\ — > ujtt and a\ — > pn can be a sensitive 
discriminator among models. They found that if the qq are created in the 3 Po configuration, as commonly assumed 
in flux-tube models, the ratio of D/S ratios is 



-(ax -» pw) _ 1 
-> ujtt) ~ 2' 



L4 = _2, (14) 

an ju . . \ 

as 



They found the same ratio in the case of qq creation by gluon exchange in the static limit ( "colour coulomb" ) but 
that it departs from —1/2 in the case of transverse gluon exchange. It was suggested that this might be useful as a 
signature of the one-gluon exchange component in the physical decay amplitude, and noted that experimentally the 
ratio is —0.35 ± 0.090, 2a away from —1/2. Today the ratio 

22. (ai -► cm) 

^4- -4 = -0.39 ±0.06 (15) 

as v 1 1 

has a greater precision (20| due to recent experiments [2l| . though the statistical deviation from —1/2 remains similar. 
Although the authors of ref [§] speculated that the common ratio for 3 P and coulomb-gluon cases is because of a 
lack of spin- flip, which is violated in the case of transverse gluon exchange and hence the deviation from —1/2 in that 
case, they did not explicitly demonstrate the source. 
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In the factorisation scheme, the amplitude for these decays is proportional to a single matrix element; this is an 
example of case 2 cited above. Once again, zeroes in the Wigner coefficients reduce the expansion of eqn. ([5]) to a 
simpler form and enforce the conservation of angular momentum (Z12 = and If = L), whereby the amplitude in a 
partial wave L is proportional to a unique spatial matrix element Al- The two decay modes of interest differ in the 
spin quantum number of the initial state, s = (61) and s = 1 (a\). The amplitudes are given by 



a.CPr ^ Sl %) = ^(-) i+1 n s { \ I I } { f f x ) 2 \ 



This gives 



as 



1 



P 1 ^ 3 S 1 1 S ) = -^=A S a D ( 1 P 1 ^ 3 S 1 1 S Q ) = -^=A D (17) 
a s ( 3 Pi ~> 3 Sj %) = -±A S a D ( 3 P! ^ 3 Si %) = ^=^L D (18) 



Thus we have established that the ratio eqn. (|14p is an immediate result of the factorization and qq creation with 
spin 1. A deviation from this ratio is indicative of a breaking of factorization, such as by a transverse gluon which 
transfers spin and momentum ( "spin-orbit coupling" ) in general. We shall return to this mechanism for breaking of 
factorization in section III. 



C. 5* and D wave decays of vector mesons 

If both of the vector states are 3 Si then the decay amplitudes are V + (0 + , 1 + , 2 + ), and the amplitudes in S, D- 
wave are each proportional to a unique spatial matrix element As, Ad] this is an example of case 3 discussed earlier. 
Decays of vector mesons provide a range of tests of factorization and decay dynamics. Decays of the type 

3 Si -> 3 Po,i,2+ 3 S! (19) 

3 Si 3 Pi.2+ 1 S (20) 

3 Si -> 1 Pi+ 3 Si (21) 

3 S X -> ^i+^o (22) 

all belong to the special case 3 described earlier; their decay amplitudes in a partial wave L are each proportional to 
a unique matrix element Al. Substituting into eqn. ([5]) I — 0. s — 7 ' — 1 for the initial state and l\ = 1 and I2 = for 
the final states gives the amplitude % 12 l for the decay in a partial wave L with final states coupled to ji2- 



^'12, l( Si ► 1 Pj 1 + S2) — y ) ni SlS2Sl2Sl2 j- l3 - I2 

S12 



(«. 


S12 




/ S12 


1 


3X2 \ 


1 312 


jl 




I L 


1 


1 J 




The results are shown in Table [J below. The pattern of amplitudes is realized in specific model calculations that 
have implicitly assumed factorization, e.g.[l3l|. which give explicit expressions for the spatial dependences Agr n \ and 
Ao(n) f° r radial excitations n. The amplitudes in Table [J differ from those in Ref. 13( by a phase associated with 
the ordering of the angular momentum coupings. 



1. ?/)(n 3 Si) — > flavoured mesons 

The results of Table U can be applied immediately to -0("- 3 Si) — > Dq^D* and also to D\D^*' . In the latter case 
data may be used to determine the mixing angle between 1 Pi and 3 Pi. 

The eigenstates for axial flavoured mesons are in general mixtures of the 3 Pi and 1 Pi states. Ref. [2^ defines the 
mixing angles by 

\D 1L ) = cos0| 1 P 1 )+ sin^Px) 

\D 1H ) = -sm^Pi) +cos0| 3 Pi) (24) 

and discusses ways of determining them experimentally. The amplitudes for axial meson production as a function 
of mixing angle follow from Table U with careful treatment of phase conventions for the spin-mixed states. Ref. [23[ 
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3 Si ^ 3 Po 3 Si 


ais = 


-As/2 




aiD — 





3 Si -^ 3 Pi 3 Si 


ais = 


A s /V3 




am = 


A D /4V3 




aw = 


Ad/A 



3 Si ^ 3 P 2 3 Si 


ais = 







am = 


A D /AVE 




0.2D = 


-Ad/aVs 




asD = 


-Aosjl/Vh 




0.3G = 





3 Si Pi 3 Si 


ais = 


A s /V6 




am = 


-A D /2V& 




CL2D — 


-A d /2V2 


3 Si -+ 1 Pi % 


a\s = 


As/273 




am = 


A D /2VS 



3 Si ^ 3 Pi % ais = A s /V& 

am = -A S /2V6 

3 Si ^ 3 P 2 'So a 2 D = A D /2V2 

TABLE I: Decay amplitudes a jl2 L for the decays (fl9|) - ([2"2 



adopted the following conventions: for qc states (as opposed to cq states) and with orbital and spin angular momentum 
combined in the order (I <X> s)j, the heavy quark limit gives <j> = —54.7° [23| so that the states are 

./2 |3 



\Dil) = \I^?i)-\J^?i) 
2,i d , , ./l,3 T 



= VS 1 1>+ V 3 1 (25) 

The amplitudes of eqn. ([5]), shown in Table HI are for the topology in which the created q (q) ends up in the meson 
with quantum numbers siliji (s 2 /2j2)- If the axial states are labelled with the quantum numbers SxZiJi they are qc 
states in accordance with conventions of ref [23| . However, the conventions in the present paper are that amplitudes 
apply to meson spin coupling in the order (s g) l)j, so there is relative minus sign associated with the 3 Pi part of the 
amplitude. Thus for the mixed Dm, D\l states in the heavy quark limit, the amplitudes for 3 Si — > D±i,D, D\hD are 

a J12i ( 3 S! D 1L D) = y|a iiai ( 3 Si Pi %) + \J\a ji3L ?$ x ^ 3 Pi %), (26) 

a, 12L ( 3 Si -> D 1H D) = y|a J12i ( 3 Si Pi %) - ^ 12 l( 3 Si ^ 3 Pi %), (27) 

and likewise for 3 Si — > D\lD*,Di H D*. This gives the relative decay widths (up to phase space corrections) shown 
in Table HJ below. 

Hence in the heavy quark limit 

r(^(n 3 Si) -> D*D 1L ) = 2r(V>(n 3 Si) DD 1L ) (28) 
T(V(n 3 Si) -> D*Di H ) = 2r(^(n 3 Si) -> D£> 1Jf ) (29) 

as well as 

r(^(" 3 Si) ^D*D lL ) = 2r(V>(n 3 Si) -^D*D ) (30) 

r(^(A)-i?*AH) = ir(^(" 3 Si) -> £>*r> 2 ) (3i) 

r^^S^DDm) = r(V(n 3 Si) -> DD a ) (32) 

In addition there is a selection rule that the D$ is produced only in helicity or 1; i.e denoting helicity states by 
0, (±), (±±) then 
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S 2 


D 2 


D D* 


1 





D 1L D* 


2 





D 1H D* 





1 


D 2 D* 





2 


D 1L D 


1 





DxhD 





i 

2 


D 2 D 





1 

2 



TABLE II: Relative widths 3 Si — > D*-Do,i.2 ° r DDo,i,2', the states Dil,h are light and heavy axial mesons in the heavy quark 
limit. 



a(^(n 3 S 1 )^D 2 (±±)D*( T )) = (33) 

This will be derived in the next section. 

2. Helicity selection rule 

In the factorization scheme, the decay of a transversely polarised 3 Si — > 3 P2+ 3 Si, with the tensor meson maximally 
polarised along the decay axis, is predicted to vanish: 

a( 3 S 1 (+)^ 3 P 2 (++) 3 S 1 (-))=0 (34) 

This selection rule is a test of factorization; a significant non-zero strength for this helicity amplitude in a decay 
1~~ — > 1~~2 ++ signals either a breakdown of factorization or the presence of 3 D! in or of 3 F 2 in 2 ++ . The 
origin of the selection rule is most transparent if we consider the helicity amplitude structure directly. Its generality 
can then be assessed by transforming to partial wave amplitudes. 

First consider the helicity picture. The decay is 

9i94 -> [qiqz] + [<73<?4] (35) 

through the creation of 5293 (Fig[l|. Denoting fermions with S z = ±1/2 by u, d respectively, the initial 3 Si(+) has its 
qq spins oriented U1U4. The final 3 P2(++) 3 Si(— ) then has to be [M1M2] + [^3^4] with the [U1U2] also having L z = +1, 
so spin-flip is required for a non-vanishing amplitude. 



<i4 




FIG. 1: Strong decay topology. 

Spin conservation on spectator lines following the steps above, or the diagrammatic techniques of ref. 0], enable 
relations among helicity amplitudes to be calculated in such factorizing models. 

In order to expose the more general dynamics underpinning this selection rule, and to exhibit the relations among 
the various helicity amplitudes, it is useful to transform between helicity and partial wave amplitudes. As before, 
consider a state with spin j decaying to two particles with spins respectively j\ , j 2 ■ The final state can be characterised 
by quantum numbers (ji2, L) or by helicity quantum numbers (Ai, A2); the translation between the two bases, for an 
initial state with spin projection m along some axis, is given by 



|jm;Ai,A 2 ) = ^ J 2 ^ (ji2A|jiAi; j 2 - A 2 }0'A|ji 2 A; L0)\jm; j 12 L) (36) 
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This enables helicity amplitudes a\ 1 \ 2 to be written as linear combinations of partial wave amplitudes a,j 12 L- We are 
interested here in decays of the type 

V(\)^ X jA^)+V(-\2) (37) 
with Ai — A2 = A, the relation between helicity and partial wave amplitudes follows from (|3€>|) above with j = J2 = 1, 

a Xl x 2 = y|(lA|iiA i; l-A a )ois + ^/|]T<lA|. n2 A; 20)<iiaA|iiAi; l-X 2 )a jiaD + V3(1A|3A; 40><3A|iiA i; l-A 2 )a 3G 

jl2 

(38) 

The resulting relations are shown in the first column of Table IIII1 These relations apply generically to the decay of 
any vector meson to any scalar, axial or tensor meson \h ■ For decays of the type 

3 S 1 (A)^ 3 P, 1 (A 1 )+ 3 S 1 (-A 2 ), (39) 

where each of the vectors are explicitly in a 3 Si state and the Xji is a 3 Pj'i state, the amplitude is obtained by 
substituting for Oj 12 l( 3 Si -^ 3 Pj 1 3 Si) from Table HI the results are shown in the second column of Table Hill The 
selection rule (|34l) is explicit in the last line of Table IIIII and follows immediately substituting 

ais = a 3G = 0; a 1D = A d /aV5; a 2D = -A D /AV3; a 3D = -A D ^7/15; (40) 











3 Si(A) ^ 3 P il (Ai)+ 3 Si (-A3) 


h 


= 


ao,o 


= ais/VS — a 1D \j2/3 


= -A S /2V3 






"0,+ 


= ais/VS + ai D /V6 


= -A s /2-/3 


h 


= 1 


ao,o 


= 


= 






a+,o 


= ais/V6 + a 1D /\/V2 — o, 2 d/2 


= A S /3V2- A D /12 






"0,+ 


= —ais/V6- ai D /^/vl — a 2D /2 


= ~A S /3V2- A D /6 






a+- 


= ais/V^ — a 1D /V3 


= A s /3V2 - A D l\2 


h 


= 2 


0-0,0 


= -ois-v/2/15 + 2ai D /v / 15 + 303D/V35 - 2a 3G V 3 / 35 


= -Ao/2V3~ 






"+,0 


= -ais/VTO - a 1D /V20 - a 2D /y/V2 + 4a 3 n/V^5 + 2a 3G /\ 


/35 = -Ad/4, 






ao,+ 


= a ls /V30 + a 1D /V60 + a 2D /2 + 2a SD /V3E + a ZG ^3/35 


= -Ad/2^/3 






a+- 


= ais/VW— a-i D /\/Z + a 3D ^3/35 - 2a 3G / v / 35 


= -Ad/4 






a++- 


= ais/Vl + am/VW - a 2D /V& + a SD i/2/105 + a 3G /V70 


= 



TABLE III: Column 1 expresses helicity amplitudes a^ 1 x 2 in terms of partial wave amplitudes dj 12 L for decays (|37[1 of generic 
vector states. Column 2 expresses helicity amplitudes in terms of spatial amplitudes Al for decays (|39p with explicit s and I 
quantum numbers. 

The amplitude a 3 G = in 3 Po or 3 Si models since (Z12 = h = 1) ® {L = 3) can couple to (If = 2,3,4), which 
cannot couple to the I = initial state by the vector decay operators V or f . The appearance of this zero can be 
tested by measurement of the various helicity amplitudes which satisfy the linear relation 

a ,o - 2/ s/Sa+fi + 2/ Via + - = a 0: + + 1/ \/6a ++ ._ (41) 

The selection rule V(+) — * T(++)V(—) = can be violated by failure of factorization, such as when single gluon 
exchange produces the q2<?3 and flips-spin such as u\ — ► d%, or if there are 3 Di admixtures in the wavefunctions of the 
produced or initiating vector mesons. The general property that breaks factorization and mixes 3 Di components in the 
produced vector meson is essentially the same: in models the latter is generated by spin-orbit coupling, such as from 
gluon exchange [24]] . To the extent that vector mesons and e + e~ annihilation are dominated by 3 Si configurations, and 
the strong decay amplitude factorizes, the selection rule will apply. For the 0(4415), which is consistent with being 
4 3 Si [l^l, the decays ■0(4415) — » DD% have been observed by initial state radiation [2(|; our selection rule may be 
testable on the high mass side of the 0(4415) in its decays to the low mass tails of D*(2010) and -D 2 (2460) respectively. 
It can also be tested in the e + e~ continuum immediately around 4.5GeV as 3 Di contamination is expected to be 
minimal 12 511 . 
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D. Hybrid and Exotic Charmonium 



While 3 Si ip states are expected to dominate the couplings to e + e~ annihilation, there are local cc resonance 
structures in in the 4-5GeV energy range whose structure is still unestablished [251] . In particular there are the enigmatic 
structures y(4260) and Y (4325) [27}. These have no natural assignment within cc spectroscopy and explanations 
include hybrid charmonium, or molecules (e.g. by either cqcq tctraquarks or DD\ and D*Dq attractive forces via tt 
exchange), or even effects associated with S-wave charmed meson thresholds [28| . 

These states are near to the S-wave thresholds for DDi, -D*-Do.i,2- Such decays are an example of case 1 discussed 
earlier: each S-wave amplitude is proportional to a single spatial matrix element. The coefficient is a function of the 
spin and orbital angular momentum of the vector initial state, and thus the pattern among decay amplitudes differs 
for 3 Si, 3 Di and hybrid interpretations, where for the latter the cc have I = 1 and s = 0. Eqn. (J5j) gives the coefficient 
of the spatial matrix element, and the results arc shown in Table [TV] For axial mesons the amplitudes are shown in 
both the 1 Pi- 3 Pi basis and in the heavy quark limit, where for the latter the amplitudes are given by cqns. (|26p and 
(|27[) and their analogues. 







3 Sl 


3 Dl 


inpl 


D D* 




-1/2 





1/3V2 




D* 


1/V6 


-1/2VE 


1/2^/3 


£i( 3 Pi 


D* 


i/Vs 


1/4^/3 


1/2V6 


D 2 D* 









1 fi 
6 y 2 


tfiCPi 


D 


l/2>/3 


l/2y/3 







D 


i/Vg 


-1/2VE 


-1/2V3 


D 1L D 




1/2 





-1/3V2 


D 1H D 







1/2V2 


1/6 



TABLE IV: Relative S wave amplitudes for vector charmonia decays with 3 Si, 3 Di and 1 IIPi (hybrid) configurations; the states 
Dil, Dih refer to axial mesons in the heavy quark limit. 



If production of charmed mesons in the decays of 3 Si cc is confirmed to factorisc, then using Table ITVl the relative 
decay amplitudes to DDi, D*-Do,i.2 may be used to determine the structure of cc states that are near to the S-wave 
thresholds. In particular this applies to Y"(4260) and F(4325). There are characteristic zeroes that may occur for 
vector meson decays: 

r( 3 Si -> D 1H D) = (42) 

r( 3 D! -> D 1L D) = (43) 

r( 1 IIPi (hybrid)] — > Di( 1 Pi)D) = (44) 

The first pair of zeroes arise from the affinity of light and heavy Dil, Dih for S and D couplings respectively, and 



the zero (|4"2"|) was noted by ref. 25j. For the hybrid decay the result follows from the conclusion of lattice QCD, 
section IIA, that decays are driven by qq creation in spin-triplet, which implies that a pair of spin-singlets (such as 
D and 1 Pi) cannot be produced from a spin-singlet, such as a hybrid vector cc . In practice these predictions will be 
affected by mixing, which can be determined from other processes (e.g. see j22l | ), a nd by phase space. The relative 
rates are insensitive to form factor effects at low momenta (see for example refs [3, [TH, |22| ) . 



III. ELECTRON-POSITRON ANNIHILATION 



We consider now the production of meson pairs in e + e annihilation, supposing that such processes proceed through 
the strong decay of a virtual quarkonia state 

e + e~ -> qiq 4 -> q{q 2 + <73<7 4 - ( 45 ) 

An analagous model for ip decays to li ght flavour meson pairs was found to be consistent with data assuming the q{q i 
state is some radial excitation n 3 Si [20|. If the same applies here, the relative production amplitudes of x Pi 1 So, 
: Pi 3 Si, 3 Pj 1 So and 3 P 3 3 Si will have the pattern of Table I, independently of n. Such relations apply in the 
limit of equal momentum decays and provided there is not an unfortunate double conspiracy in which both a single 
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n dominates and a node in its amplitude coincides with the kinematic region of interest. This may be checked by 
varying q 2 to see if the ratios are stable or vary in an oscillatory or nodal manner. The assumption that the pair 
q{q~4 is dominantly in a 3 Si configuration is reasonable above charm threshold where the coupling e + e~ — » 3 T>\(cc ) 
is theoretically and empirically suppressed [2fl |25| . The results of Table U then apply immediately to e + e~ — ► D*Dq ^ 
and also to D^Di, which in the latter case may be used to determine the mixing angle between 1 Pi and 3 Pi. While 
this is strictly true on a 3 Si ip resonance, it may also be expected to hold through the 4-5 GcV region of interest 
where 3 Si is predicted to dominate the e + e~ cross section. 

Application to e + e~ — > ip + x j a l so follows if this is dominated by strong flux-tube formation and breaking. We 
shall argue in section IIIB that this is more likely to be dominated by (perturbative) gluon exchange, which breaks 
factorization and gives a different pattern of amplitudes than strong flux-tube breaking. Our results may be used to 
test this hypothesis. 

In the case of light flavours the neglect of e + e~ — > 3 Di(qq ) is more questionable. The leptonic widths of 3 Di(qq ) 
are nonetheless expected to be relatively small (30j. and empirically the known vector mesons app ear to fit well with 
(radially excited) 3 Si with some mixing with hybrid vectors without need for significant 3 Di[31|. This is clearly an 
area whose phenomenology merits further clarification. To that end we apply our results with 3 Si dominance to 
light flavours in the hope of shedding further light on this sector and isolating 3 Di states. For ip decays this simple 
assumption appears to be consistent with existing data[20l. [32| . 

A. e + e~ — > flavourless mesons 

In the case of e + e~ — ^neutral states, charge conjugation restricts the production of axial- vector mesons in association 
with h or 1" to e + e~ — > V + 3 Pi or 1 So + 1 Pi . The amplitudes of Table [J apply and the relative rates then follow 
by application of equation ([5]) : 

3 Si+ 3 P : 3 Si+ 3 Pi : 3 Si+ 3 P 2 : % + 1 Pi = 3S 2 : AS 2 + D 2 : 6D 2 : S 2 + D 2 (46) 

and hence 

a( 3 Si + 3 Pi) = ^ CT ( 3 Si + 3 P ) + ia( 3 Si + 3 P 2 ) (47) 
3 6 

together with 

3a( 1 S + 1 Pi) - <7( 3 Sr + 3 P ) + ia( 3 Si + 3 Pa) (48) 

and their corrolary 

aCS Q + 1 Pi) = ia( 3 Si + 3 P 2 ) + l ( x( 3 S 1 + 3 Pi). (49) 

Note that necessarily 

a{e+e- ^ 3 Si + 3 Pi) > a(e + e~ ^ 3 Si + 3 P ). (50) 

For flavoured states the two axial mesons are mixtures of *Pi and 3 Pi; whatever the mixing angle may be, the 
inequality holds true in the sense that the 3 Si + 3 Po production rate cannot exceed those of both of axial mesons. 
In the case of charge conjugation eigenstates we are restricted to applying it to light flavours or to e + e~ — > ?/; + Xj- 
The former case is less well controlled theoretically, due to relativistic effects and potential contamination from 3 Di 
background in e + e~ annihilation, though the above relations appear to be consistent with data and are discussed 
in ref. f32|. One of the central applications of the present paper will be to test these predictions against data on 
e + e~ — ► ip + Xj where preliminary indications are that the relation eq. (|50p is violated (33|. This is discussed in section 
IIIB. 

The amplitude V(—)T(++), where here T denotes a tensor meson, should also be measured for light flavours 
where VT modes are prominent, especially in ip decay. Within the factorization hypothesis and 3 Si dominance the 
V f<x cannot be produced with / 2 maximally polarised; a[e + e~ — > V(j z — —l)f%{jz = +2)] = . This may be 
studied in e + e _ — > 5ir — 27r + 27r _ 7r° by isolating the channel the ui being a narrow state can enable the angular 
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distribution in decay — > tt + tt~ to be measured. The main background here is the potential contamination from 
e + e~ annihilation in the 3 D! state. Although models and data do not suggest this is significant, nonetheless one 
cannot rule it out. If the amplitude is empirically found to be small, in accord with the selection rule, one could turn 
this to advantage and study the amplitude as a function of q 2 and observe if it turns on in the neighbourhood of the 
predicted 3 Di resonances around 2.2GeV[30l]. 

B. Factorization breakdown and preformation by OgE 

Data on e + e~ — > ip + X at 10.6 GeV cm. energy show three prominent enhancements X in e + e~ — -> ip + X 33], 
which are consistent with being the r} c , r]' c and xo- The observed pattern of states appears radically different to what 
is seen for light flavours, for example the apparent prominence of e + e~ — > ip + xo with only a hint of xi an d much 
suppressed X2 contrasts with light flavours where e + e~ — ► ufc is clearly seen[20|. This suggests that this process for 
heavy flavours may be controlled by a production mechanism where factorization is broken. 

On theoretical grounds one expects that strong factorization may be overwritten here. In e + e~ annihilation at 
E > 6 GeV, creation of an initial cc leaves up to 3 GeV available. As the cc separate, forming a strong flux tube 
up to 0(lfm) long, the energy of O(lGeV) enables light-flavoured qq to form. That is the familiar dynamics that 
appears to be realized at low energies for light flavours. In the present example, the most probable circumstance is 
that the excess energy produces multiple qq leading to final states DDttit.... The experimental selection on final 
states ipX(cc) isolates an unlikely configuration where the 3 GeV has produced a cc exclusively. For the flux tube 
to grow without splitting until it contains 3 GeV of energy would require it to extend to distances exceeding A.qc D ■ 
This is exponentially unlikely with increasing energy. 

Alternatively the energy can be transmitted through a single gluon which converts to cc . While this is pcrturbative 
and expected to be sub-dominant for processes involving light flavour creation, the question arises at what energy or 
for what flavours this dominates over flux-tube breaking. The purpose of this section is to propose ways of answering 
this by experiment. We make specific reference to e + e~ — * (cc) + (cc) as there are emerging data in the form of 
e+e - -> V + X. 

As momentum flows through the gluon, it can transfer spin or angular momentum between the cc to which it is 
coupled. In general therefore we anticipate that factorization will break down. 

Ref Q have considered these matrix elements in the explicit non-relativistic limit - (see Appendix B of ref 0], 
especially eqs B5-B7). In that limit the gluon-exchange operation transforms as S.S and L.S but there is no S.L 
operator (where the first operator refers to the transformation property of the gluon emission and the second operator 
to that of qq creation). Thus in the strict non-relativistic limit of that model, the V(—)T(++) selection rule would 
appear to survive for the decay of a 3 Si vector meson. This is no surprise following the discussion after ea.(|34p : 
non-zero amplitude requires spin flip at the emission vertex and orbital flip at the cc creation vertex; while the former 
occurs in the non-relativistic limit, the latter does not. 

However, in e + e _ annihilation at q 2 = E 2 m , the production of a cc allows an S.L operator at Q(q 2 /m 2 ) . An 
explicit calculation of the gluon exchange contributions to e + e~ — > ip + Xj nas been made in NRQCD in ref |34j | and a 
non- vanishing amplitude for V( — )T(++) is found even at threshold, in accord with the discussion above. Threshold 
is when q 2 = 16m 2 ; the amplitudes depend upon r 2 = 16m 2 jq 2 . At high energies, where r 2 — > the contribution from 
e + e~ -^ 3 Di — * cc will become increasingly important while for the threshold region, r 2 — ► 1, the e + e~ —> 3 Si — > 
cc becomes more dominant. 

At the 10.6 GeV cm. energy of the data [33j|, r 2 = 0.28, and ref [34| finds for the one-gluon exchange (OgE) 
contribution to the cross sections a(ipxo '■ V'Xi : 4'X2) ~ 12 : 2 : 3, which contradicts eq. (|5"0|) based upon factorization 
and assumption of a 3 Si initial state. In the threshold limit r 2 —> 1 the analysis simplifies and comparison between the 
predictions of gluon exchan ge a nd factorization becomes sharpest. In this limit the VT amplitudes for transversely 



polarised initial state of ref. [34j satisfy 



a[V(-)T(++)] : a[V(0)T(+)] : a[V{+)T(0)} - 1 : l/y/2 : 1/V6 (51) 

in accord with S-wave dominance and the results of Tabl e IIIII The relative cross-sections from the OgE mechanism 
for e + e~ — > ipXo,i,2 m vicinity of threshold r 2 — > 1 in ref [34l | become 

a(^Xo ■ <PX1 ■ M - 24 : 2 : 3. (52) 

Compared to the results at higher energy, r 2 = 0.28, the relative sizes of ipxi ■ "0X2 have not changed much but there 
is a significant relative enhancement of cr(ipxo) near threshold. 
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This prediction, that the cross-section for ipxo dominates, contrasts with the results of factorization near threshold. 
For 3 Si initial state in the S'-wave region near threshold 



afe) -> (53) 
3 

r 



3 

ct(V'Xo) = to-(V'Xi) 



Analogously, for a 3 Di initial state 



<7(VXo) -> (54) 
5 

c-(V'Xi) = 3 Cr (V'X2) 

which is also utterly unlike the OgE predictions. Finally one may allow for a coherent mixture of 3 Si and 3 Di initial 
state. Results become model dependent but cr(^xo) cannot be made larger than both a{ipxi) an( i cr (' t PX2)- Thus in 
the region of threshold there appear to be marked differences in the expectations of factorization, eqs (p>3")) . ([54")) and 
OgE cq.flSU). 

As one increases energy above threshold, for 3 Si initial state, the D-wave decay enables 0(^X2) to turn on but 
with the amplitude a[V(-)T(++)\ = or at least small compared to a[V(0)T(+)] and a[V(+)T(0)]. This also 
contrasts with the predictions from OgE where the [V(— )T(++)] amplitude is the largest for VT production, cq (|5Tj) . 
A possible contamination comes from e + e~ -^ 3 Di — > cc contributions which may not be negligible at 6-7GcV c.o.m 
energies. The S-wave decay amplitudes from initial 3 Si, 3 Di and also from hybrid vector mesons are compared in 
Table IIVI Above threshold where Z?-wave decays are important and 3 Si- 3 Di mixing is allowed, results are highly 
model dependent. While it may be possible to force a(ipXo) to dominate by suitable choice of mixing angle, this is 
not expected to hold true as a function of q 2 . 

Thus if dominance of ipXo is confirmed over a range of q 2 away from threshold, this would support OgE as the 
dominant decay mechanism. Conversely, if data near threshold confirm a[V( — )T(-\ — h)] — > 0, this would signal factor- 
ization being dominant. In any event, we anticipate that the relative populations and helicity structures of tpXj wn l 
vary with q 2 . We recommend that this be investigated in e + e~ annihilation at super-B factories by means of ISR to 
access a range of energies. In particular experiment should attempt to measure the spin dependence of e + e~ — ► ipX2 
as a function of q 2 and compare with the analogous amplitudes in e + e~ — > w/2. 



IV. CONCLUSION 



The factorization property of strong decay triggered by qq creation in spin-triplet, as revealed by lattice QCD, 
merits further testing. This general feature leads to relations among amplitudes, which can be used as further tests 
of this dynamics and to determine the nature of participating mesons. Thus we have identified the following tests. 

1. 0(n 3 Si) decays or e + e~ annihilation in the 4-5 GeV energy range will not produce D* D2 with the tensor meson 
in helicity two. This tests whether the dynamics revealed by lattice QCD for light mesons applies more generally 
for the strong creation of light flavours. 

2. If confirmed, then the production e + e _ — > D\ may be used to determine the axial meson mixing angles in 
the 3 Pi- 1 Pi bases. 

3. If the mixing angles are known from elsewhere, the pattern of charm pair production can identify the nature of 
the decaying ip state. This has an application of immediate relevance in determining the nature of the enigmatic 
charmonium-like structures Y(4260) and Y(4325) and also of -0(4415). Determining whether the cc content 
of these states is S = (as for a hybrid) or S = 1 then follows from the relative production rates of various 
combinations of charmed mesons, in particular of their DD\ branching ratios. 

4. The application to light flavours in e + e~ is less solid, but measurement of the w/2 amplitudes as a function of 
q 2 may isolate 3 Di resonances in the e + e~ channel. 

5. For the creation of heavy flavours, as in e + e~ — > ipXj, empirical and theoretical arguments suggest that produc- 
tion is dominated by a single hard gluon rather than the factorization mechanism. The apparent excess of ipxo 
and absence of 1/1x2 needs establishing. We expect that the pattern of Xj states and their helicity amplitudes 
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will vary significantly with q 2 . We identify the threshold region e + e~ — > ipXj between 6.5 and 7.5 GeV as 
particularly promising for determining the relative importance of single hard gluon and strong factorization for 
heavy flavours. 
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